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I review in this work the properties of nuclear and quark matter equation of state 
under the influence of a strong magnetic field. We analyse the results when they 
are applied to compact stars description. Finally, we describe the emission of 
gravitational wave from a compact star that goes over a phase transition. 

1 Introduction 

The determination of the nuclear matter equation of state and the appearence 
of new phases at high densities is one of the main goals in present nuclear 
physics. Nowadays, huge recently planned experiments have brought excel- 
lent expectations on these determinations. RHIC and LHC are the most 
promissing tools for direct observations of the high density matter behaviour, 
from heavy ion collisions. Another important experiments arc linked to the 
observation of neutron star properties, once they are believed to behave like 
giant hypcrnucleus and/or giant quark bags. In this direction, Chandra X- 
Ray Observatory and BeppoSax have done an excellent job analysing the 
data originated by X-Ray or thermal emission from pulsars. We emphasize 
the data that indicates important results about more massive neutron stars 
(M~2.0-2.2M©) and others about smaller stars with radius ~ 30% smaller 
then ordinary neutron star radius. New astrophysical tools are related to 
the emission of gravitational waves from these extremely dense objects and 
shall be observed by the american LIGO or by the french-italian collaboration 
VIRGO. 

Such a variety of experiments being located for the seek of nuclear mat- 
ter equation of state and the quark-gluon plasma phase transition implies in 
a large ammount of theoretical works which are capable to improve the ex- 
periments capabilities, direct them to better targets and have their payback 
with results that may prove theories. Nowadays one can find different works 
described in many papers, not only because of theoretical physicist's cria- 
tivity but also because of the many different approachs high energy physics 
support. One can go from the complex analytical quantum chromodynamics 
(QCD), step on the lattice simulations, drive through feynman diagrams, and 
get wherever he wants. They all are validy and the scope of this work do not 



GFMarranghello: Workshop SMFNS/ICIMAF - February 2, 2008 1 



include a discussion of the best way of doing nuclear physics or the best model 
to fit the bulk properties of matter. We intend to use a quantum hadrody- 
namics (QHD P based model and the MlM bag model to describe compact 
stars main properties and trace a consistent path from bulk nuclear matter 
description to neutron star properties and the gravitational wave emission. 

2 Nuclear Matter Equation of State 

One consistent, well based and effective model to describe nuclear matter is 
based on the pioneered work by Walecka which consists one of the basis of 
nuclear matter effective field theory. Further developments were done by many 
authors from whose we extract the works by Boguta and BodmeJ™ Zimanyi 
and MoszkowskPl and Taurines et. al0, developed in the following years. 
We have used the model as described by the lagrangian density below 
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which includes the baryonic octet (p, n, A,S + , S°, S~, 3°) coupled to 
mesons (<j,u>, g), and free lcptons (e, /i). The scalar and vector coupling con- 
stants, g,j, go, and the coeficients b, c are determined through the nuclear 
matter bulk properties fitting, i. e., binding energy E& (= -16.3 MeV), com- 
pression modulus K (= 240 MeV) and effective nucleon mass M* — M — g a a 
(= 732 MeV) at saturation density po (=0.153 fm~ 3 ). The g e coupling con- 
stant is determined through the nuclear matter asymme try c oefficient , a4 (= 
32.5 MeV). All the coefficient values are described in ref!^liJ 

Applying standard technics of field theory and solving the equations at a 
mean-field level one can extract expressions for the thermodynamical quanti- 
ties of pressure, p, energy density, e etc. and verify the suitability of this model 
results to the experimental values of compression modulus, K and effective 
nucleon mass, M* , resulting in excellent results for neutron star masses and 
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radii, which can be compared to previous work J ' '. In figure ^ we describe 
the behaviour of the density of each different particle described by the model 
as a function of the total baryon density. 




P, (&>') 

Figure 1. Particle distribution as function of baryon density. 



3 Phase Transition to Quark Matter 

The nuclear to quark matter phase transition is expect to occur at high den- 
sities and/or high temperatures, restoring the chiral symmetry. We follow 
in this section the work by Heiselberg and Hjorth-JenserP whose perfectly 
described the conditions for the occurence, or not, of a mixed phase com- 
posed by hadrons and deconfined quarks simmultaneously. First claimed by 
GlendenninglQI the occurence of a mixed phase where quark structures are 
immersed into a hadron matter (or vice- versa at higher densities), is now 
uncertain. This work will not clarify this topic, once it is model-dependent, 
however, as an application of this study we intend to propose the possible 
difference on the dctcctation of such phases inside hybrid stars. 

Generally, the properties of the transitions with one or more than one 
conserved charges are quite different. The most important feature is that 
pressure may be constant or vary contineously with the proportion of phases 
in equilibrium. Reviewing Heiselbergs work, one can search the answer of 
which kind of transition is more energetically favored through an analysis of 
the bulk energy gained in such a transition and the Coulomb and the surface 
energy relation, as we stress below. 

Assuming the model described in the previous section and the MIT bag 
model for the quark matter, 
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where f2 is the thermodynamical potential for the MIT bag model, we relate 
the energy densities in order to calculate the bulk energy. Giving the large 
uncertainties in the estimates of bulk and surface properties, one cannot claim 
that droplet phase is favored or not, once it depends crucially on nuclear and 
quark matter properties. 

If droplet sizes and separations are small compared with Debye screening 
lenghts, the electron density will be uniform to good approximation, and the 
Debye screening lenght Ac is defined by 



(3) 



where n,, in and Qi are the number density, chemical potential and charge of 
particle species i. 

In the system we consider, screening effects can be estimated and, if the 
characteristic spatial scales of structures are less than about 10 fm for the 
nuclear phase, and less than about 5 fm for the quark phase, screening effects 
will be unimportant, and the electron density will be essentially uniform. In 
the opposite, the total charge densities in bulk nuclear and quark matter will 
both vanish. 

When quark matter occupies only a small fraction, / = VqmVvnm ' °^ 
the total volume, quarks will form spherical droplets. The surface energy 
per droplet is given by eg = aAirR 2 , where a is the surface tension, and the 
Coulomb energy is 



ec = -|^-(/0qm - Pnm) 2 R 5 - (4) 

Minimizing the energy density with respect to R one obtain the usual result 
e s = 2ec and find a droplet radius 
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Our results, considering the models described before, cannot be definitive 
before a high precise determination of the surface tension, a. However, for 
a ~ 70— lOOMeV, a mixed phase is unfavored and the star will have a density 
discontinuty. The result is a larger radius to support the mixed phase and 
the equation of state for the constant pressure phase transition is presented 
below. 



E (fin ) 

Figure 2. Equation of state for pure neutron star matter (solid line) and for the transition 
to quark matter (dotted line). 



4 Magnetic Field Effects 

Several indenpendent arguments link the class of Soft j-Ray Repeaters (SGRs) 
and perhaps anomalous X-ray pulsars with neutron stars having ultra strong 
magnetic fields. In addition, two of four known SGRs directly imply, from 
their periods and spin-down rates, surface fields in the range 2 — 8 x 10 14 G. 
The magnitude of the magnetic field strength needed to dramatically affect 
neutron star structure and an estimate done by Lai and Shapiro yields 
B~2x 10 18 ^/^rjtmf a m ^ ne m t er i° r °f neutron stars. 

In this section we investigate the effects of very strong magnetic fields 
upon the equation of state of dense matter in which hyperons and quarks are 
present. In the presence of a magnetic field, the equation of state above nuclear 
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saturation density is significantly affected both by Landau quantization and 
magnetic moment interactions, but only for fields strengths B > 1 18 G. 

The previously presented lagrangian density in section|21can be re- written 
in order to include the magnetic field as 



B 
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For detaile d and well founded texts on magnetic field effects the authors 
recommend ref! ^^ , from where one can find the energy spectrum for the 
protons given by 
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In a strong magnetic field, particles motion is perpendicular to the field 
lines and are quantized to give discrete Landau orbitals, and the particles 
behave like a one-dimensional gas rather than a three-dimensional gas and, 
the higher is the field, the lower is the number of occupied Landau levels The 
energy of a charged particle changes significantly in the quantum limit if the 
magnetic field H > H c . For electrons, ff c ~ 4 x 10 13 G, u and d (massless 
quarks), H c ~ 4 x I0 15 G, s-quark, H c ~ I0 18 G and protons, 10 20 G. 

Considering the phase transition aspects, a strong magnetic field unfavor 
the mixed hadron-quark phase once it increases the eletron number in the 
nuclear phase leaving the quark phase praticaly unchanged. This i ncr eases 
the droplet size as it was discussed in the previous session and in ref.^. 
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5 Compact Stars 



The first step through an analysis of compact stars recall the simpler solution 
to the general relativistic Einsteins equations which represents the static and 
spherical symmetric stars, the Schwarschild solution, written as the Tolman- 
Oppenheimer-Volkoff equation d 1 ^^ ! The TOV equations describe the struc- 
ture of a static, spherical and isotropic star wit the pressure p(r) and the en- 
ergy density e(r) reflecting the underlying nuclear model. The TOV equations 
involve various constraints and boundary conditions: they must be evaluated 
for the initial condition e(0) = e c (with e c denoting the central density) and 
M(0) = at r — 0; the radius R of the star is determined under the condition 
that on its surface the pressure vanishes, p(r)\ r= R = 0. 

The condition for chemical equilibrium for neutron stars are: 

Hi = biH n - qi(pi) (7) 

where \ii and fii stand for the baryon and lepton chemical potentials, respec- 
tively; bi is the baryon number and the baryon and lepton electrical charges 
are represented by 

The corresponding equations for baryon number and electric charge con- 
servation are: 



Pbary 



Elf. <■» 



3tt 



and 



B I 

Here one can visualize the importance and the effects of quark and nu- 
clear matter parameters determination, once they will have reflects on the 
properties of compact stars. First, quark stars maximum mass and radius are 
directly governed by the value of the bag constant as: 

^ _ 1.964M ^IQgfrn (1Q) 

V -^60 V -^60 

where B§o = B / (6QMeV/ fm 3 ) in the massless quarks case. In the same sense, 
the strong coupling constant a c is related to the quark star properties. Both 
B and a c will determine whether or not not the hadron-quark phase transition 
will take place. The sequence of neutron stars (mass-energy density relation) 
and the energy density in the interior of a neutron star is presented in the 
figures below. 
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Figure 3. Masses of neutron star as a function of central energy density. The hybrid star 
without mixed phase (solid line), with mixed phase (dotted line) and the pure neutron star 
(dashed line) are presented. 
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Figure 4. Energy density distribution inside the neutron star for two different central en- 
ergies. 



The nuclear matter coupling constants gsa, 9Buj and gs e present the 
same sensibility and uncertainties. Nuclear matter results for effective mass 
compression modulus and particle distribution play also an important role on 
the phase transition and star properties and, consequently, on the emission of 
gravitational wave as we will discuss in the following sections. 

The magnetic field effects are presented in the following figures of equation 
of state, neutron star masses and particle distribution. 
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Figure 5. Equation of state at B=0 (solid line) and B=10 18 G with (dotted line) and without 
(dashed line) the effects of anomalous magnetic moment. 
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Figure 6. Masses of pure neutron star at B=0 (solid line) and B=10 18 G with (dotted line) 
and without (dashed line) the effects of anomalous magnetic moment and the masses of 
hybrid stars at B=10 18 G (long dashed line) and B=5 X 10 18 G (dot-dashed line) as functions 
of central energy density. 



6 The Conversion of Neutron Stars to Hybrid Stars 

The transition from configuration H to HQ may occur through the formation 
of a mestastable core, built up by an increasing central density. The increase in 
the central density may be a consequence of a continuous spin-down or other 
different mechanisms the star could suffer. This transition releases energy, 
exciting mainly the radial modes of the star. These modes do not emit GWs, 
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Figure 7. Particles distribution at B=0 (solid lines) and B=10 18 G (dotted line) as functions 
of the baryonic density. 



unless when coupled with rotation, a situation which will be assumed here. 

Strange matter is assumed to be absolutely stable and a seed of strange 
matter in a neutron star would convert the star into a hybrid or strange 
star. The speed at which this conversion occurs was calculated by O untold, 
taking into account the rate at which the down- and strange-quark Fermi seas 
equilibrate via weak interactions and the diffusion of strange-quarks towards 
the conversion front. 

Accordingly to^-^, cold neutron stars can convert to strange star matter 
with speeds ranging from 5km/ s to 2 x I0 4 km/s. The outcome of such an 
event would emit an incredible amount of energy in 0.5ms to 2s. We apply 
this discussion to the partial conversion of a neutron star to strange matter, 
forming a hybrid star. 

However, the emission of a gravitational wave, which will be discussed 
in the next section, can only recognize the phase transition after a structural 
rearrangement of the star, which shall reduce its radius and gravitational 
mass, conserving its tot al baryon number. Such mechanism has already been 
studied by Bombaci for the emission of 7-ray bursts. 



7 Gravitational Wave Emission 

In order to simplify our analysis, we will consider that most of the mechanical 
energy is in the fundamental model. In this case, the gravitational strain 
amplitude can be written as 

h(t) = hoe-W^- 1 " ^ (11) 
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where ho is the initial amplitude, u>q is the angular frequency of the mode and 
r gw is the corresponding damping timescale. The initial amplitude is related 
to the total energy E s dissipated under the form of GWs by the relational 

1/2 

(12) 



, 4 
h = 

w r 



GE g 



where G is the gravitational constant, c is the velocity of light and r is the 
distance to the source. 

Relativistic calculations of radial oscillations of neutron star with a quark 
core were recently performed by Sahu et al!^l However, the relativistic mod- 
els computed by those authors do not have a surface of discontinuity where 
an energy jump occurs. Instead a mixing region was considered, w here the 
charges (electric and baryonic) are conserved globally but not locally^l Os- 
cillations of star models including an ab rup t transition between the mantle 
and the core were considered by Haensef^l MiniuttP^ and Marrang hellcP. 
However a Newtonian treatment was adopted and the equation of state used 
in the calculations does not correspond to any specific nuclear interaction 
model. In spite of these simplifications, these hybrid models suggest that 
rapid phase transitions, as that resulting from the formation of a pion con- 
densate, proceed at the rate of strong interactions and affect substantially 
the mode frequencies. However, the situation is quite different for slow phase 
transitions (the present case), where the mode frequencies are quite similar 
to those of a one-phase star^! In this case, scaling the results of 1 ^ 1 , the 
frequency of the fundamental mode (uncorrected for gravitational redshift) is 
given approximately by 

~(M/M e ^ 1/2 



i/ « 63.8 



R 3 



kHz (13) 



where the mass is given in solar units and the radius in km. 

Once the transition to quark-gluon matter occurs, the weak interaction 
processes for the quarks u, d and s 

u + s — ► d + u and d + u — ► u + s (14) 

will take place. Since these reactions are relatively slow, they are not balanced 
whi le t he oscillations last and thus, they dissipate mechanical energy into 
heaM According to calculations of reference [21], the dissipation timescale 
can be estimated by the relation 

nn Vl50MeF\ 4 (M Q \ 

^°- oi (^H s (15) 
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where m s is the mass of the s-quark in MeV and M c is the mass of the 
deconfined core in solar masses. This equation is valid for temperatures in 
the range 10 8 — 10 9 K. On the other hand, the damping timescale by GW 
emission is 

where again the stellar mass is in solar units, the radius is in km and the 
rotation period P is in milliseconds. 

In a first approximation, the fraction of the mechanical energy which will 
be dissipated under the form of GWs is f g = ^ 1+T X ■ Notice that the 
damping timescale by GW emission depends strongly on the rotation period. 
Therefore one should expect that slow rotators will dissipate mostly of the 
mechanical energy into heat. In table 4 is given for each star model the 
expected frequency of the fundamental mode (corrected for the gravitational 
redshift) , the critical rotation period (in ms) for having f g = 0.50, the GW 
damping for this critical period and the quality factor of the oscillation, Q 

—TTV T gw . 

Considering now a strongly mangetized neutron star whose B ~ 10 18 G, 
we detect an increasing in the gravitational mass of a star composed by pure 
hadronic matter. However, the gravitational mass of hybrid stars decrease 
for the case of increasing magnetic field and the explanation is quite simple: 
once the hybrid star has lower gravitaional mass than the pure hadronic one, 
due to the presence of a softer quark matter in its inner shells and, the phase 
transition from hadron to quark matter occurs at lower densities with higher 
magnetic fields, the highly magnetized hybrid star will have a more important 
quark-gluon plasma core and a lower mass. This fact will represent a larger 
energy released during the phase transition and gravitational wave emission 
which could be detected from larger distances. 

After filtering the signal, the expected signal-to-noise ratio is 



DC 



(s/Ny = 4i 4tt^ ( 17 ) 

JO &n(y) 

where h(v) is the Fourier transform of the signal and S„(^) is the noise power 
spectrum of the detector. Performing the required calculations, the S/N ratio 
can be written as 

\2 _ ^ 1,2 ( T gw \ Q 2 



{S/N) ' = l hi {s^)Tfw (18) 
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Table 1. Oscillation parameters for B=0: the damping timescale r gw is given for the critical 
period; maximum distances for VIRGO (V) and LIGO II (L) are in Mpc, and the LIGO II 
result for B=10 18 G. 





p 

1 CTlt 


T gw 


Q 


^max 


^rnax 




(kHz) 


(ms) 


(ms) 




(VIRGO) 


(LIGO II) 


B = 10 18 G 


1.62 


1.64 


87.0 


442 


4.9 


10.2 


11.5 


1.83 


1.25 


27.0 


155 


6.4 


13.5 


15.28 


2.06 


1.13 


17.0 


110 


6.0 


12.8 


14.48 


2.32 


1.06 


11.5 


84 


5.1 


11.1 


12.56 


2.72 


1.00 


8.4 


72 


3.6 


5.7 


6.45 



In the equation above, the angle average on the beam factors of the detector 
were already performed. 

From eqs.( 4) and (11), once the energy and the S/N ratio are fixed, 
one can estimate the maximum distance T> max to the source probed by the 
detector. In the last two columns of table 4 are given distances D ma x derived 
for a signal-to-noise ratio S/N = 2.0 and the sensitivity curve of the laser 
beam interferometers VIRGO and LIGO II. In both cases, it was assumed 
that neutron stars underwent the transition having a rotation period equal to 
the critical value. 

We emphasize again that our calculations are based on the assu mption 
that the deconfinement tr ans ition occurs in a dynamical timescale \ZQ In the 
scenario developed in refffSl a mixed quark-hadron phase appears and the 
complete deconfinement of the core occurs according to a sequence of quasi- 
equilibrium states. The star contracts slowly, decreasing its inertia moment 
and increasing the its angular velo city until the final state be reached in a 
timescale of the order of 10 5 yr^HI Clearly, in this scenario no gravitational 
waves will be emitted and this could be a possibility to discriminate both 
evolutionary paths. 

8 Additional Effects 

Additional effects and/or scenarios may improve, or not, the detection of 
gravitational waves from compact stars. An important aspect is related to 
the colour superconducting phase. As, in principle, this newer phase would 
turn the equation of state even softer, enlarging the energy release during the 
astrophysical phase transition, we decide to analyse such a phase. 
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The superconducting quark phase is described by the thermodynamical 
potential 



n C FL = Sl F ree - — o A 2 /X 2 + B . (19) 

The resulting equation of state was evaluated in the TOV equations lead- 
ing to a smaller radius star, as expected, but also to a less massive star, 
even for the gravitational and for the baryonic mass. This results did not 
change significantly the properties of gravitational wave emission as we have 
expected, however, we still expect a better understanding about high-density 
matter properties to better describe such phase. 

The convertion from a neutron star to a strange star, instead of a hy- 
brid star, is a more catastrophic event and may release much more energy, 
being detec tab le from longer distances. This kind of transition was studied 
by Bombaci^, in order to describe the emission of 7-ray bursts. Applying 
our results to describe the emission of gravitational waves we got the fol- 
lowing results which leads to one of the most promissing gravitational wave 
generators. 

Table 2. Properties of a NS— >SS conversion and the maximum distances for VIRGO (V) 
and LIGO II (L) are in Mpc 



M bar 


M ss 


Rss 


y\ V irgo 
max 


T\Ligo 
max 


1.0749 


0.6445 


8.40 


15.47 


34.63 


1.3202 


0.7804 


8.87 


20.01 


44.79 


1.5724 


0.9241 


9.30 


24.27 


54.33 


1.8342 


1.0727 


9.68 


28.20 


63.13 


2.1045 


1.2210 


9.98 


32.07 


71.79 



9 Conclusions 

Summarizing, we worked on a nuclear many-body theory with parametrized 
coupling constants and the MIT bag model. The phase transition from hadron 
to quark matter was described considering screening effects and the equation 
of state was determined. We applied the resulting EOS to the TOV equations 
also considering intense magnetic field effects. These results are according to 
previous results obtained by many authors. As a tool for gravitational wave 
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emission, we considered the phase transition occuring inside the compact star 
obtained from the TOV equations and releasing energy. According to our 
results, such an event would be detected as far as 13 Mpc without considering 
magnetic fields and 15 Mpc for compact star with B=10 18 G. Considering a 
NS— >SS transition, we obtained a much larger distance of up to 71 Mpc. 

References 

1. B.D. Serot, J.D. Walecka, Adv. Nucl. Phys. 16:1,327,1986; Int. J. Mod. Phys. 
E6:515-631, 1997. 

2. A. Chodos, R. L. Jaffe, K. Johnson, C. B. Thorn and V. F. Weisskopf, Phys. 
Rev. D9, 3471 (1974). 

3. J. Boguta and A. Bodmer, Nucl. Phys. A292, 413 (1977). 

4. J. Zimanyi and S. A. Moszkowski, Phys. Rev. C42, 1416 (1990). 

5. A.R. Taurines, C.A.Z. Vasconcellos, M. Malheiro, M. Chiapparini, Phys. Rev. 
C63:065801, 2001. 

6. G.F. Marranghello, C.A.Z. Vasconcellos, J. A. de Freitas Pacheco, M. Dillig, 
Int.J.Mod.Phys. Ell:83-104, 2002. 

7. G. F. Marranghello, C. A. Z. Vasconcellos and J. A. de Freitas Pacheco, Phys. 
Rev. D66, 064027 (2002). 

8. N. K. Glendenning, F. Weber and S. A. Moszkowski, Phys. Rev. C45, 844 
(1992). 

9. H. Heiselberg and M. Hjorth- Jensen, nucl-th/9902033 

10. N. K. Glendenning, Phys. Rev. D46, 1274 (1992). 

11. D. Lai and S. Shapiro, ApJ. 442, 259 (1995). 

12. A. Broderick, M. Prakash and J. M. Lattimer, astro-ph/0001537 

13. A. Broderick, M. Prakash and J. M. Lattimer, Phys. Lett. B, 531 (2002) 167.. 

14. D. Bandyopadhyay, S. Chakrabarty and S. Pal, astro-ph/9703066 

15. R. C. Tolman, Phys. Rev. 55, 364 (1939). 

16. J. R. Oppenheimer and G. M. Volkoff, Phys. Rev. 55, 374 (1939). 

17. V. Olinto, Phys. Lett. B, 192 (1987) 71. 

18. I. Bombaci, Lect. Notes in Phys. 578:Physics of Neutron Star Interiors, edited 
by D. Blaschke, N. K. Glendenning, A. Sedrakian (Springer- Verlag, 2001). 

19. J. A. de Freitas Pacheco, Proceedings of the VI International Workshop on 
Relativistic Aspects of Nuclear Physics, ed. T.Kodama et al., World Scientific, 
2001, p.158 

20. P.K. Sahu and G.F. Burgio and M. Baldo, ApJ 566, L89 (2002). 

21. P. Haensel and J.L. Zdunik and R. Schaeffer, A&A 217, 137 (1989). 

22. G. Miniutti and J. A. Pons and E. Berti and L. Gualtieri and V.Ferrari, 
astro-ph/0206142 

23. Qing-de Wang and Tan Lu, Chin.Astr.Astrophys. 9, 159 (1984). 

24. P. Haensel and M. Prozynski, ApJ 258, 306 (1982). 

25. G. Lugones and J. E. Horvath, Phys. Rev. 66, 074017 (2002). 



GF Marranghello: Workshop SMFNS/ICIMAF - February 2, 2008 15 



